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Abstract 

The spectral functions of the vector current and the axial-vector current have been mea¬ 
sured in hadronic r decays using the OPAL detector at LEP. Within the framework of 
the Operator Product Expansion a simultaneous determination of the strong coupling 
constant a s , the non-perturbative operators of dimension 6 and 8 and of the gluon con¬ 
densate has been performed. Different perturbative descriptions have been compared 
to the data. The Contour Improved Fixed Order Perturbation Theory gives a s (m 2 ) = 
0.348±0.009 ex p±0.019 t heo at the r-mass scale and a s (m|) = 0.1219±0.0010 exp ±0.0017theo 
at the Z°-mass scale. The values obtained for a s (m|) using Fixed Order Perturbation 
Theory or Renormalon Chain Resummation are 2.3% and 4.1% smaller, respectively. 
The ‘running’ of the strong coupling between so — 1.3 GeV 2 and so = m 2 has been tested 
from direct fits to the integrated differential hadronic decay rate R t (sq). A test of the 
saturation of QCD sum rules at the r-mass scale has been performed. 
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1 Introduction 


The r lepton is the only lepton heavy enough to decay into hadrons. A comparison of the 
inclusive hadronic decay rate of the r with QCD predictions can give fundamental parameters 
of the theory. The energy regime governed by m T = 1.777 GeV is regarded as a compromise 
region between the low and high energy regimes where non-perturbative and perturbative QCD 
dominate, respectively. In fact, r decay is probably the lowest-energy process from which the 
coupling constant a s can be cleanly extracted W without large complications from non- 
perturbative effects, while the perturbative expansion still converges well. 

In this analysis the most important quantity to measure is the strong coupling constant 
a s {m 2 ). The ‘running’ of a s , for energy scales smaller than m T , can be tested with the integrated 
differential decay rate into hadrons di? T /ds, where y/s denotes the mass of the final-state 
hadronic system and R r = T(r —> h v T )/V{r —> ev e v T ) is the hadronic decay width of the tau 
normalized to the decay width of the tau going into electron and neutrinos. This is possible as 
the hadronic decay rate R T (s 0 ) depends on the strong coupling constant a s (s 0 ) at the scale s 0 
only, where So denotes the upper integration limit for the integral over di? T /ds. 

The measured a s {m 2 T ) can be transformed into a value for a s (m|) through the renormaliza¬ 
tion group equation (/3-function). In doing that, the relative error of a s (s) decreases like the 
decrease of a s (s) itself. After the evolution to the Z° mass the strong coupling is reduced to 
a s (m|) ~ (l/3)a s (m^) and its error is reduced to Aa s (m|) ~ (l/9)Aa s (m^). Hence, the signif¬ 
icance of this measurement compares favorably with other a s (m|) determination methods ||. 

Inclusive observables like the hadronic decay rate R r (s o) have been calculated in perturba¬ 
tive QCD to O(ctg). Some remaining theoretical uncertainties due to corrections in powers of 
1/ml can be avoided if the differential decay rate di? r /ds is measured and compared to the 
theory by means of its spectral moments which are weighted integrals over d_R r /ds. As a result, 
the power corrections and a s can be simultaneously determined from a fit. While R r [m 2 ) can 
be precisely determined from the leptonic branching ratios and the r lifetime, dR r /ds involves 
a measurement of the invariant mass of the hadronic system. Thus, an exclusive reconstruction 
of all hadronic final states in r decays is necessary. 

In this paper an analysis is presented using data taken with the OPAL detector at LEP at 
energies within ±3 GeV of the Z° peak. The analysis includes measurements of the differential 
decay rates diiQv/A/ds for vector (V) and axial-vector (A) decays and their respective spec¬ 
tral moments. Using these moments, fits of QCD predictions are made extracting the strong 
coupling constant a s (m 2 ) and parameters of the non-perturbative expansion, most notably the 
contributions of dimension 6 and 8 operators. The measurement is based on a set of spectral 
moments defined by the same weighting functions used by ALEPH [[7], |J and CLEO 0 . 

The differential decay rates themselves can be re-expressed in terms of spectral functions 
of the vector and axial-vector current, r;(s) and a(s). This measurement serves for saturation 
tests of QCD sum rules at the r-rnass scale by comparing the experimental values of the sum 
rules with chiral QCD predictions. Furthermore, by evaluating the moment integrals between 
zero and So, where is an energy smaller than m T , the ‘running’ of a s is tested in a single 
experiment. 

The theoretical framework for inclusive observables from hadronic r decays is described 
in section ^]. After a short description of the OPAL detector in section [| the selection of 
hadronic r decays is described in section f|. In section |] the unfolding procedure is described. 
The measured and unfolded spectra are discussed in section [| followed by a description of the 
systematic uncertainties in section |7]. Section |8| contains the results for the moments of R r and 
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for the spectral functions. The extraction of the strong coupling constant and of the power 
corrections, from fits to the moments of R r (s 0 ), is discussed in sections |9| and |T0|. Section 0 
describes the test of the ‘running’ of ct s . The application of QCD sum rules to the spectral 
functions is discussed in section O. Finally, the results are summarized in section O. 


2 Theoretical description of hadronic r decays 

QCD predictions of inclusive observables in hadronic r decays have been calculated including 
perturbative and non-perturbative contributions. These observables can be related to the differ¬ 
ential, non-strange hadronic decay width, normalized to the decay width of t~ —> e~v e v T 



2 r 


1 + 2 ^, 


Imi 7 V/A( S ) 


+ Imi 7 V/A( S ) 


( 1 ) 


where s denotes the square of the invariant mass of the hadronic system and the labels V 
and A stand for the vector and axial-vector contributions, respectively^. Sew — 1-0194 is an 
electroweak correction term |T(| and |l/ ut i| 2 = 0.9512 ± 0.0008 is the squared CKM weak mixing 
matrix element |IT[. The functions ImlJ are proportional to the spectral functions for the 


non-strange currents with angular momenta J — 1 and J = 0 as indicated by the superscripts. 
The latter spectral function vanishes for the vector current, since no scalar particle has been 
observed in r decays, while Irn/7 A is given by the pion pole, assuming that the pion is the only 
pseudo-scalar final-state in non-strange r decays: 


ImJ7»( s ) = 



e^e^r) Ay 


( 2 ) 


with N 1 r being the number of selected r decays into pions. The spectral functions for the vector 
and the axial-vector currents are defined in equation (|22|). 

Within the framework of QCD weighted integrals or moments of ([I]) have been calcu¬ 

li 

S ° ' ' U (_S_\ l dRrWA 

S 0 / V rn r / 


lated 


rjfcZ 

n r,V/A 


(•So) = J ds 


l-± 


ds 


( 3 ) 


The moments are used to compare the experiment with theory. In what follows, ten moments 
for kl = 00,10,11,12,13 for V and A are used. The first moments i?(?° v ^ A ('m 2 ) are the total 
normalized decay rates of the r into vector and axial-vector mesons given by ([I]) integrated over 
s. In the naive parton model these two rates are identical and add up to the number of colors. 
Since only non-strange currents are considered in this work the naive expectation has to be 
multiplied by |Kid| 2 - Including the perturbative and non-perturbative contributions, equation 
(H) is usually written as ||T2| : 


<v/ a(»o) = |s EW | V„d| 2 I 1 + O»o) + C.N) + E 

Z ' D= 2 , 4 , 6 ,... 


r D,kl/ \ 

"v/a( s o) 


( 4 ) 


where S'ew is the same multiplicative correction as in equation (HD and are additive elec¬ 
troweak corrections. The latter has been calculated for kl = 00 only |D| yielding <5 ew ( m l) — 


1 The notation V/A will be used throughout the paper to indicate vector and axial-vector contributions, 
respectively. 
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_ q. 0010. In the higher moments it is assumed that this term scales with the integral 
over the weight functions in equations ([!]) and (j3j) like the 0(a s ) correction: 


X ,kl 

°EW 


( s o) — 


^r!naive( S o) 


sr/00 


R 00 .. 

''t nni 


(Til , 
u r, naive V" v tJ 


2\^Ewl m: ? 


( 5 ) 


Therefore the contribution to the moments is small (~ 0.1 %) and the uncertainty due to 
this term is neglected in the analysis. The other factors in equation (j4|) are explained in more 
detail below. 


2.1 Perturbative correction terms d kl , 

pert 

The perturbative term ^ rt is known to third order in a s 
in a s [O. For kl = 00 and So = m 2 it is: 


C.K) = ^ 

7r 


h5.2023 Q ' s +26.366(78-003 +K 4 ) ^^- +0 (a 5 s (m 2 r )) . (6) 

7T" 7 T 6 7T 4 V 7 


,« 3 K) 


and partly known to fourth order 




This result which truncates after the fourth power of a s is refered to as Fixed Order Perturbation 
Theory (FOPT). Different attempts have been made to obtain a resummation of some of the 
higher order terms. The resummation scheme proposed in |F|] compensates for higher order 
logarithmic terms in a s by expressing the <5pe rt (so) terms by contour-integrals in the complex 
s-plane along the circle |s| = So and solving numerically for each a s (s) along the circle (Contour- 
Improved Perturbation Theory, CIPT). The different a s (s) values on the circle can be calculated 
from a s (m 2 ) by solving numerically the /3-function: 


da 

dins 


= / 3(a) = —ftia 2 — /3 2 a 3 — /? 3 a 4 — // 4 a 5 + 0(a 6 ), 


( 7 ) 


with a = a s (s)/7r, /3i = 9/4, f3 2 = 4, P^ s = 10.0599 and = 47.2306 [[L4J for 3 quark flavors. 
The last two coefficients are renormalization scheme dependent and the quoted values belong 
to the MS-scheme. The third method considered in this paper resums the leading term of the 
//-function to all orders in a s by inserting so-called Renormalon Chains (RCPT) 


One of the leading theoretical uncertainties for FOPT and CIPT comes from the unknown 
O (ctf) correction Ji 4 . Expanding the perturbative corrections in terms of CIPT gives: 


l+Ct( S ») = E^n(«o). 


( 8 ) 


n>0 


where the functions A f p are the weighted contour integrals. For kl = 00 the function is: 




ds ( a s (—s) 


f R 4 A 4 N 

n pi pv p^ 

2 4 - 2 4 + 4 - 2 ± + 24 - 4 


7T 


mi 


mi 


mi 


mi 


mi 


mi 


( 9 ) 


In the MS-scheme and for three flavors the first four terms are: K 0 = Ki = 1, K 2 = 1.63982, 
K 3 = 6.37101 Hl8l-fm A bold guess for A' 4 gives A' 4 ~ K 3 (K 3 /K 2 ) ~ 25 |12|. Similar estimates 
are given in |23|. [21|| . A central value of K 4 = 25 is used, with an uncertainty of AA' 4 = ±50 in 
the perturbative expansions for CIPT and FOPT. 


2 The fixed-order corrected version (up to the third order in a s ) quoted in the lower portion of table 6 in 

ref. M is used in the fit. 
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Another major theoretical uncertainty is the choice of renormalization scale fi in the a s 
dependence of £pe rt (so)- The scale ratio ( = fi 2 /s 0 is varied from 0.4 to 2.0 in all three models 
described above as suggested in [|J. 

The choice of the renormalization scheme (RS) can also alter the result. Following the 
prescription in || the third coefficient of the /3-function 0. j/ s is varied between 0.0 and 2.0 /3^ s 
in order to obtain the uncertainty due to different renormalization schemes. 


2.2 


Power correction terms 


In the framework of the Operator Product Expansion (OPE) f25[ the non-perturbative contri¬ 
butions are expressed as a power series in terms of 1/m 2 absorbing the long-distance dynamics 
into vacuum matrix elements (0(0)) |2fil [28 


Thus, they can be written as sums over power 


corrections of different dimensions, D: 

s:kl 

"non-perpV/A 


( s o) - X) 

D= 2,4,6,.. 


<-D,kl 

°V/A 


( s o)- 


( 10 ) 


In contrast to the perturbative part described in the previous section the power corrections 
differ for the vector and the axial-vector currents. 

In equation © the correction of dimension D = 2 is a mass correction term and therefore 
belongs to the perturbative part. The D = 4 term is the Erst term with major non-perturbative 
contributions, namely the quark condensates for the three light flavors (0ip) u ,d,s and the gluon 
condensate GG). If one neglects the small s-dependence of the power corrections, the <5yp 
terms can be expressed for all kl values by a product of the same (vector/axial-vector) operator 
of dimension D (or the power correction for kl = 00) and a simple integral over the ^/-dependent 
weight-functions |T!| : 


r D,kl 
°V/A 


(so) 


D=2 

D =4 

D =6 

D=S 

D =10 


kl 

( 1 

0 

-3 

—2 

0 

\ 

00 

1 

ml 

so 

-3 

—2 - 3— 

So 

—2 — 


10 

0 

-1 

m T 

3 

2 + 3— 


11 

0 


so 

o 

SO 



0 

l 


-3 


12 


0 

0 

so 

-l 

so 

) 

13 


x y- £v/a(aO ( 0 ( 0 )) 


( 11 ) 


where each entry in the matrix belongs to a particular dimension D and a particular moment kl, 
as denoted by the first row and the last column. The parameter jj, is an arbitrary factorization 
scale which separates the long-distance non-perturbative effects, which are absorbed in the vac¬ 
uum matrix elements (0(0)), from short-distance perturbative effects which are incorporated 
in the Wilson coefficients Cy/ a(0) |E^ . 


This approach is used for the dimension D = 6 and D = 8 terms, taking Sy^°° as free 
parameters. For the dimension D = 2 and D = 4 terms the full s-dependence is taken into 
account for the theoretical description of the moments [12]]. The least precisely known D = 4 
parameter, the gluon condensate, which is known only to 50% j4|, is also taken as a free 
parameter in the fit, while the D = 2 term is calculated from the quark masses and the strong 
coupling. 

Terms with dimensions higher then 8 are neglected in this analysis as they do not contribute 
to R00 /a as can be seen from equation (pT|) . 
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3 OPAL detector 


29|. A brief description of the 


A detailed description of the OPAL detector can be found in 
features relevant for this analysis follows. 

A high-precision silicon microvertex detector surrounds the beam pipe. It covers the angular 
region of | cos0| < 0.8 and provides tracking information in the r-<p (and 0 after 1992) direc¬ 
tion^ |3^,|3l|. Charged particles are tracked in a central detector enclosed inside a solenoid that 
provides a uniform axial magnetic field of 0.435 T. The central detector consists of three drift 
chambers: a high-resolution vertex detector, the large-volume jet chamber and the ^-chambers. 
The jet chamber records the momentum and energy loss of charged particles over 98 % of the 
solid angle and the ^-chambers are used to improve the track position measurement in the z 
direction [Q . 

Outside the solenoid coil are scintillation counters which measure the time-of-flight from 
the interaction region and aid in the rejection of cosmic events. Next is the electromagnetic 
calorimeter (ECAL) that is divided into a barrel (|cos0| < 0.82) and two endcap (0.81 < 

| cos 6*| < 0.98) sections. The barrel section is composed of 9440 lead-glass blocks pointing 
to the interaction region. Each block subtends approximately 10 x 10 cm 2 with a depth of 
24.6 radiation lengths. The two endcap sections consist of dome-shaped arrays, each having 
1132 lead-glass blocks, mounted coaxial with the beam, where each block covers 9.2 x 9.2 cm 2 
with a typical depth of 22 radiation lengths. The hadron calorimeter (HCAL) is beyond the 
electromagnetic calorimeter and instrumented with layers of limited streamer tubes in the iron 
of the solenoid magnet return yoke. In the region | cos 6 1 1 < 0.81 this detector typically has 
a depth of 8 interaction lengths. The hadron calorimeter is covered by the muon chamber 
system, composed of four layers of drift chambers in the barrel region and four layers of limited 
streamer tubes in the endcap region. 


4 Event selection and reconstruction of r decays 

OPAL data collected from 1990 to 1995 is used in this analysis. The data were taken within 
±3 GeV of the Z°resonance. The Monte Carlo samples used in this analysis consist of 600 000 
r-pair events generated at y/s = mz with Koralz 4.0 [|33j. Their decays were modelled with 


Tauola 2.4 and then processed through the Geant [|35| OPAL detector simulation |36 
The non-r background Monte Carlo samples consist of 1 000 000 qq events generated with 
Jetset 7.4 P7fl , 800 000 Bhabha events generated with Radbab 2.0 QQ], 600 000 fi -pair 


events generated with Koralz 4.0 33] and 800 000 events from two-photon processes generated 


with Vermaseren 1.01 40, 41 . 


4.1 Selection of r-lepton candidates 


The standard r selection procedure as described in |42|] begins with the rejection of cosmic 
rays, multi-hadronic events and events from two-photon processes. Cosmic rays are rejected 
by the time-of-flight information of the tracks. Multihadrons are removed from the sample by 
requiring two narrow jets (cones with a half opening angle of 35 °) and up to six tracks in the 
event. The events from two-photon processes are eliminated by allowing an acollinearity angle 
of up to 15 ° between the two jets. 

3 In the OPAL coordinate system the cc-axis is horizontal and points to the center of LEP. The y-axis is 
vertical and the 2-axis is in the e _ beam direction. The angle 9 is defined relative to the 2-axis. 












































The remaining event sample contains tau pairs, Bhabha events, and muon pairs. Events 
with an energy deposit of more than 0.8 x 2E^, eam are identified as Bhabhas. An event is 
classified as a muon pair if two tracks carry energy of more than 0.6 x 2E' beaxa and if both 
tracks have at least two hits in the muon chambers and almost no energy deposit in the ECAL. 
The remaining events are classified as r pairs if the polar angle of the total cone momentum 
calculated from track momenta and ECAL clusters satisfies | cos#| < 0.95 for both cones. 

After this selection both cones in each event are treated independently. The non-tau back¬ 
ground is further reduced by requiring one or three tracks in each cone with a total charge 
of plus or minus one. A total of 297 988 r candidates survive these selection criteria with an 
estimated non-r background fraction of 3.9 %. 


4.2 Identification of r-decay modes 


A Maximum Likelihood selection as used in previous publications (see e.g. j43|) is applied 
to the data and the Monte Carlo samples to distinguish between the following decay modes: 
t~ — > z/ T X“, where X - is one of e~77 e , /x - F M , n~, 7i~7i 0 , 7i~2n 0 1 7r~37r°, 27r“7r + , 27r _ 7r + 7r°, 
27r _ 7r + 27r°. The charge and parity conjugated modes are implicitly assumed for r + —> F T X + 
decays. Fourteen reference distributions are used to distinguish between the different one-prong 
channels and five reference distributions are used in the three-prong case. Decays with charged 
kaons instead of pions are suppressed by a cut on the specific energy loss dE/dx in the drift 
chamber. Decays into electrons are distinguished from the other modes by the ratio E/p of the 
ECAL energy associated with the cone over the track momentum, and the dE/dx information. 
Muons are identified by the number of hits in the muon chambers and the outermost HCAL 
layers. The different hadronic decay modes with zero, one, two or three neutral pions are 
separated by using the number of reconstructed photons in the ECAL (see section fL3|) . 

The decay channels used in this analysis are the three non-strange one-prong modes with 
at least one neutral pion: Tin 0 , 7r27r° and 7r37r°, and the three non-strange three-prong modes: 
37r, 37T7T 0 , 37r27r°. A total of 65899 r candidates are selected^ in these channels with an esti¬ 
mated background fraction of 26.6 % including misidentihed r decays and the remaining non-r 
background fraction of 0.8%. Details about the treatment of the cross-talk between the signal 
channels due to misidentihed r decays are subject to sections and || 

The most important observable for the discrimination between vector and axial-vector chan¬ 
nels is the number of neutral pions in a cone. A new method to reconstruct neutral pions in r 
decays has been employed, which is described in the following section. 


4.3 Reconstruction of neutral pions 

Neutral pions are identified by their decay into two photons. Since photons are only detected 
in the ECAL, an iterative fit of photon energies and directions to the observed energy deposits 
in the ECAL blocks is performed. 

The energy deposition in an ECAL block can be expressed as a function of the photon 
energy and the photon direction. This is done by parameterizing the integrated energy density 
of an electromagnetic shower for each ECAL block. In the barrel region of the ECAL where the 
blocks have a quasi-pointing geometry only lateral shower profiles need to be parameterized. 
They can be approximated by the sum of two exponential distributions representing core and 

4 For the decay mode 7T7T 0 a cut on |cosd| < 0.9 is used in order to reduce the background from Bhabha 
events. 
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halo components [44]. For the endcap region, where the blocks are oriented parallel to the 
beam, lateral and longitudinal profiles are important. The longitudinal profile is reasonably 
described by the gamma distribution |[I] . 

The mean energy deposit of a minimum ionizing particle is subtracted from all ECAL blocks 
hit by a charged particle. The fit then Ends the smallest number of photons needed to explain 
the measured energies and provides their corresponding three-vectors. 

Energy depositions from hadronic interactions of charged pions in the ECAL are accounted 
for by assigning photon candidates which are close to track intersections with the ECAL to 
the track. The maximum angle allowed between a photon candidate and a track to which the 
photon candidate can be assigned depends on the polar angle of the track and varies between 
1.2° and 1.7° in the barrel region and between 2.0° and 3.4° in the endcap region. A photon 
candidate close to a track is still classified as a photon if the total energy of photon candidates 
assigned to this track exceeds the measured track momentum. 

All possible two-photon combinations are then used to hnd n° candidates. The combination 
resulting in the largest number of 7T° candidates with an average invariant mass deviation from 
the 7T° mass less than 1.5 a is selected. The error on the invariant two-photon mass, a, is 
calculated from the error matrices of the above photon fit. The n° four-momenta are then 
calculated from the energies and directions of the photon pairs of the selected combination 
after a constrained fit to the 7r° mass. Figure [I] shows the photon-pair mass for selected 7T7r° 
candidates before the 7r°-mass constraint. 


For all the one-prong modes a minimum energy of 0.7 GeV for each reconstructed 7T° 


is 


required while E n o > 2.0 GeV is required in the three-prong modes to suppress fake 7r°’s intro¬ 
duced by the energy deposition of charged pions in the ECAL. 



yyMass (GeV) 


Figure 1: The j^-mass in the 7T7r° channel for decays with two reconstructed photons with 
a minimal energy of 0.5 GeV. OPAL data is shown as data points; the total Monte Carlo 
prediction is given by the open histogram and the shaded histogram denotes the r and non-r 
background. 

The granularity of the ECAL allows the reconstruction of both photons from a 7T° only if 
its energy is below 12 GeV. The 7r°’s with larger energies have photons which are merged in 
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the ECAL. Therefore all the photon candidates with energies above this value are considered 
to be 7T°’s and their momentum is calculated from the reconstructed energy corrected by the 
mass of the 7r°. 


5 Unfolding procedure 


The Monte Carlo predictions for the measured spectra and their background contributions are 
corrected with the most recent constrained branching ratios of the r given in [[□]]. Effects 
due to limited detector resolution and efficiency are accounted for by a regularized unfolding 
technique 

To unfold measured distributions in s (the squared mass of the hadronic final state) the 
detector simulation is used to create response matrices which map the generated distribution 
in x = Strue to a y = s meas distribution one would measure including all detector effects. The 
following convolution integral describes the general relation between a true distribution f(x ) 
and a measured distribution g(y): 


-t-max 

g(y)= J &xA(y,x)e(x)f(x) + b(y), 


( 12 ) 


where A(y, x) is the detector response function, b(y) denotes the background distribution, and 
e(x) is the selection efficiency. Equation (|l2j) can be simplified to a matrix equation of the form: 

g = A • f + b. (13) 

This is performed in two steps. First, the true distribution f(x) is parameterized with a set of 
m parameters fj and m basis functions Pj(x) which are defined below in equation 


f(x) = J2fjf^c(x)Pj{x), 

3 =1 

with f M c{x) being the generated Monte Carlo distribution. By defining: 

^max 

Aj(y)= f dxe(x)fuc(x)A(y,x)p j (x), 


(14) 


(15) 


equation (12) takes the form: 


a(y) = ^2fjMy) + b (y)- 

3 = 1 

In the second step, g(y), b(y) and Aj(y) are represented by n bins: 

Vi Vi Vi 

<h = I Ayg(y); b,= / dy6(y); A tj = I d yAj{y). 


(16) 


(17) 


Vi — 1 Vi— 1 Vi— 1 

The basis functions Pj(x) used in equations ([14]) and (]T^) are chosen as cubic B-splines and 
thus have the following form: 


Pj( x ) = g x < 


z 3 


2 = 

(x — tj )/d 

tj ^ x tj+\ 

(1 + 3(1 + z( 1 

- Z))z) 

z = 

(x - t j+1 )/d 

tj+ 1 A x < tj. j_2 

(1 + 3(1 + z( 1 

1 

1 

z = 

(x - t j+2 )/d 

tj+2 < X < tj .|_3 

1 

CO 


z = 

(x - t j+3 )/d 

tj- 1-3 A X tj- 1_4 

0 




otherwise 
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where d = (x max — x min )/(rn — 3) is the distance between adjacent knots tk = ic min + (k — 4) d 
for k — 1,..., m — 1 knots and m splines. 

The coefficient vector f = {fj} is now observed in a fit to the data bins gi, and the unfolded 
result can be obtained with equation fll4]) . This particular choice of basis functions and normal¬ 
ization leads to the simple prediction fj = 1 for all j if the Monte Carlo generated distribution 
and the unfolded result are identical. 

In certain cases unfolding produces results with unphysical behavior. Statistically insignif¬ 
icant components of the fitted coefficient vector f can lead to large oscillations of the unfolded 
distribution. Therefore the unfolding needs to be modified by a regularization step which sup¬ 
presses the statistically insignificant parts of the solution. This is achieved by applying a smooth 
damping function to the unfolded result. The magnitude of the fluctuations is measured by the 
total curvature r(f) of the function /(x)//mc(^) : 


r(f) = 



•^min 


( d 2 f(x) 
\da: 2 fucix) 



f T ■ C ■ f, 


(19) 


where C is a constant, symmetric, positive semidefinite matrix obtained from the second deriva¬ 
tives of the basis functions pj. The regularized result is now obtained by adding the total 
curvature r(f) weighted with a regularization parameter p to the y 2 in the fit and minimizing 
the sum: 

1 

Xreg 


2 (f) = y 2 (f) + ^r(f). (20) 


The final unfolded distribution in s is given by weighting the Monte Carlo distribution in s true 
with the regularized coefficient vector f obtained from the fit (14). 

This method is bias-free as long as the detector simulation is correct for all s and is inde¬ 
pendent of the used Monte Carlo distributions provided that only the statistically insignificant 
components of the fitted contributions are damped. Possible biases due to the detector simu¬ 
lation are accounted for in the systematic errors as described in section f7|. The correct choice 
of the regularization parameter p can be tested in the following way: 

a) The coefficients fj, which are correlated and have in general different errors, can be 
transformed into a set of independent parameters a 3 which have unit variance and are 
sorted with the regularization measure pr( f) in order of decreasing significance j45|. All aj 
with j > n(p), where n(p) is the number of effective coefficients remaining after damping 
with the parameter p, have to be consistent with zero. 

b) Furthermore the y 2 -probability of the fit without the regularization term r(f) should 
increase for the regularized coefficient vector f. 


The regularization parameter is chosen according to these criteria. 

The background, in a particular one-prong (three-prong) channel, consists mainly of misiden- 
tihed other one-prong (three-prong) r decays, introducing correlations between the spectra. In 
order to provide a proper treatment of the correlations, the three one-prong (three-prong) 
channels are unfolded simultaneously. In addition to the three detector response matrices for 
the three one-prong (three-prong) signal modes, six more detector response matrices are used, 
mapping the Monte Carlo generated distribution in s true of a background channel to the back¬ 
ground part introduced by this channel in the s meas distribution of a signal channel as correlated 
background. Non-r background and other misidentihed r decays are treated as uncorrelated 
background b * as described above. 
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6 Discussion of the measured spectra 

The selection efficiencies and background fractions from simultaneously unfolded channels (cor¬ 
related background) and other r- and non-r-background sources (uncorrelated background) are 
listed in table [1|. Figures [2] and [| show the measured s meas distributions of the six channels 


channel 

efficiency 

correlated 

background 

uncorrelated 

background 

selected 

decays 

7T7T U 

(28.7 ±0.1)% 

(7.7 ±0.2)% 

(7.9 ±0.1)% 

32316 

7r27T° 

(18.8 ±0.1)% 

(45.0 ± 0.6) % 

(8.4 ±0.1)% 

13814 

7T37T 0 

(8.0 ±0.2)% 

(70.0 ± 2.2) % 

(11.4 ±0.5)% 

1738 

3n 

(34.6 ± 0.1) % 

(9.7 ±0.3)% 

(3.8 ±0.1)% 

14321 

37T7r 0 

(11.0 ±0.1)% 

(21.3 ± 1.0)% 

(6.1 ±0.3)% 

2455 

37t27t° 

(8.3 ± 0.4) % 

(82.3 ±2.9)% 

(7.1 ±0.5)% 

1255 


Table 1: Efficiencies, background fractions and total number of selected r decays. 


used in this analysis in comparison to the fitted signal after the regularized unfolding, and the 
Monte Carlo predictions. 

The 3n spectrum shows a significant deviation from the shape predicted by the Monte Carlo 
(the dashed histogram) as has been observed in previous analyses of the 3tt decay current |HB| . 
There is also a slight deviation on the left side of the peak in the 7T7r 0 channel and in the upper 
tail region. The other modes are statistically consistent with their Monte Carlo predictions. 

The x 2 values for the one-prong and three-prong fits after the regularization step are 
Xi-pr./d.o.f- = 94.0/109 and x!_ P r,/d.o.f = 71.4/69 leading to the ^-probabilities 0.85 and 
0.40, respectively. 

The unfolded distributions of the measured spectra are shown in figure f|. The plotted data 
points are strongly correlated due to the unfolding procedure. The deviations from the Monte 
Carlo prediction seen in figures |2] and j| are still present after the unfolding, most prominently 
in the 7T7t 0 and the 3 tt channel. The enhancement in the upper tail (see figures (b) and 1(a)) 
of the 7T7r° distribution can be explained within the Kiihn-Santamaria model |I7| by enlarging 
the fraction of p(1450)’s and p(1700)’s in the p decay amplitude. A similar correction to the 
three-pion current, modelled as a Breit-Wigner decay chain ai — > pn —» 37T [[IT] 

Carlo, does not account for the observed discrepancy. 


in the Monte 


7 Systematic uncertainties 

Possible origins for systematic effects on the reconstructed value for the squared hadronic 
mass, Smeas, come from the uncertainty in the energy scale for reconstructed photons and the 
uncertainty in the momentum scale for tracks, while the wrong choice of the regularization 
parameter p in the unfolding can distort the unfolded distributions. 

The energy resolution can be tested by measuring the invariant mass of the two photons 
from 7r° decays. A systematic shift in the observed mass in the data compared to the detector 
simulation can be translated into a scale factor for the reconstructed photon energies. Devia¬ 
tions of (0.5 ± 0.9) MeV for m n o have been observed between data and Monte Carlo (figure jTJ) . 
This corresponds to an energy scale factor of 1.004 ± 0.007. The energies of the reconstructed 
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• OPAL 
□ 3i jt° fit result 
Bl 31 231°, 3t 33t° bkg. 
I other bkg. 
tau MC 



s (GeV ) 



Figure 2: The measured s meas spectra for 1-prong decays. Plots (a) and (b) are the 7T7r° chan¬ 
nel, (c) and (d) are the tt2tt 0 and tt37t° modes, respectively. The points denote OPAL data 
(statistical errors only). The open histograms show the fitted spectra after the regularized un¬ 
folding, refolded into detector space. The background contributions from simultaneously unfolded 
channels (correlated background) are shown as light grey areas while the background from other 
sources (uncorrelated background) is represented in dark grey. 
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Figure 3: The measured s meas spectra for 3-prong decays. Plot (a) is the 3 tt channel, (b) and (c) 
are the 3mi 0 and 3n2ii 0 modes, respectively. The points denote OPAL data (statistical errors 
only). The open histograms show the fitted spectra after the regularized unfolding, refolded into 
detector space. The background contributions from simultaneously unfolded channels (correlated 
background) are shown as light grey areas while the background from other sources (uncorrelated 
background) is represented in dark grey. 
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Figure 4: The unfolded s true spectra. Shown are the three vector channels (left) and the three 
axial-vector channels (right) together with the Monte Carlo prediction. There are strong cor¬ 
relations between the data points due to the unfolding. The plots (a),(d),(e) are the unfolded 
spectra of plots (a),(c),(d) in figure [| and the plots (b),(c),(f) are the unfolded spectra of the 
plots (a),(b),(c) in figure [|. The error bars include statistical and systematic uncertainties. 
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photons in the Monte Carlo samples are varied by ±0.7% in order to estimate the systematic 
error due to this effect^. 

The uncertainty in the momentum of the tracks have been studied using /i pairs. The Monte 
Carlo is corrected for observed deviations between data and Monte Carlo in the mean and the 
width of the momentum distribution. The momenta and the momentum resolution of all tracks 
in the Monte Carlo samples are scaled due to the uncertainties in these corrections, thus leading 
to the quoted systematic errors. 

The damping parameter p in the regularization step of the unfolding procedure is calculated 
from the number of effectively remaining spline coefficients after the regularization. This num¬ 
ber is chosen so that the test conditions a) and b) given in section are satisfied. The default 
value (16 effective splines from 48 total splines for the 1-prong fit and 16 effective splines from 
36 total splines for the 3-prong fit) is varied by ±4 for both fits, where the range is derived from 
Monte Carlo tests of the unfolding procedure: the tests consist of unfolding fake data samples 
in the p —* m r° channel. The mass and the width of the p in the fake data samples at the 
generator level are different from the values used in the standard Monte Carlo which is used to 
create the response matrix. The allowed range of the damping parameter is then determined by 
comparing the unfolded fake data samples with their generator level distributions for different 
choices of the damping parameter, for which the test conditions a) and b) are satisfied. Within 
this range the unfolded distributions reproduce the mass spectrum of the modified p without 
biases towards the generator distribution of the standard Monte Carlo. The uncertainty due to 
the variation of the damping parameter is added as a systematic error on the unfolded results. 

Uncertainties of statistical nature from the errors on the branching ratios (see table ^]), the 
limited statistics of signal and background Monte Carlo samples, and on the efficiencies are 
incorporated in the unfolding procedure by adding them in quadrature to the statistical errors 
on the data. 

Systematic effects related to photon and 7r° detection efficiency are largely covered by the 
systematic errors. 


8 Results 


8.1 Moments of R T 

The unfolded spectra of the hadronic modes shown in figure f| are normalized to their branching 
fractions and summed up to the vector and axial-vector spectra with their appropriate weights: 



B{r —> hy /A U t ) Wy/A diVy/A 
B(t —► eu e u T ) iVy/A ds 


( 21 ) 


where iV v /a is the number of taus that decay into the hadron 1i v /a plus neutrino, and %/a 
denotes the appropriate weight of the hadronic mode to the vector or axial-vector current. 
The branching ratios of the hadronic modes (and the lepton channels), together with their 
contributing weights for the vector and axial-vector spectra, are summarized in table 

The hadronic modes urn, lothi 0 and rjTnr 0 involve decays of in’s and ff s, and do not conserve 
isospin symmetry, since their decay can occur via the electromagnetic interaction. Therefore, 
the unfolded distributions in the 37T mode, which is considered to belong to the axial-vector 

5 Since the invariant two-photon mass depends also on the angle between the two photons, this energy scale 
factor accounts for systematic uncertainties in the energy resolution and the angular resolution of the ECAL. 
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r —> t/ T X 

B 

[%] 

Wy 

w A 

comment 

eu e 

17.83 

±0.08 

- 

- 



17.35 

±0.10 

- 

- 


7T7r 0 

25.24 

±0.16 

1.0 

0.0 


37T7T 0 

4.26 

±0.09 

1.0 

0.0 

including iwtt and utitt 0 

7r37T° 

1.14 

±0.14 

1.0 

0.0 


71 

11.31 

±0.15 

0.0 

1.0 


37T 

9.26 

±0.12 

0.0 

1.0 

3h — 2Kvr — K2vr including lvtt 

7T27T 0 

9.27 

±0.14 

0.0 

1.0 


37r27r° 

0.50 

±0.05 

0.0 

1.0 

including utitt 0 and r]7nr° 

57T 

0.075 

±0.007 

0.0 

1.0 

MC 

7T47T 0 

0.12 

±0.06 

0.0 

1.0 

MC 

37T37T 0 

0.11 

±0.06 

1.0 

0.0 

MC 

57r7T° 

0.022 

±0.005 

1.0 

0.0 

MC 

KK° 

0.16 

±0.03 

1.0 

0.0 

MC 

2Ktt 

0.10 

±0.03 

0.5 ±0.5 

0.5 ±0.5 

MC 

2K°7t 

0.10 

±0.02 

0.5 ±0.5 

0.5 ±0.5 

MC 

KK°7T 0 

0.14 

±0.03 

0.5 ±0.5 

0.5 ±0.5 

MC 

LOTT 

0.21 

±0.01 

1.0 

-0.2 

MC excluding 37T7r° 

C07T7T 0 

0.046 

±0.007 

-0.25 

1.0 

excluding 37r27r° 

717171° 

0.17 

±0.03 

1.0 

-0.24 

MC 

Y 

-‘'-strange 

2.67 

±0.14 

- 

- 



Tabic 2: Branching ratios for the hadron modes and lepton channels. Shown are the fitted values 
from the Particle Data Group |71|/ and the contributing weights for the vector and axial-vector 
current. Channels marked with MC are ‘generator-level ’ Monte Carlo channels included in the 
spectra. Negative weights are used to subtract inclusively measured contributions from the wrong 
current. 


current, and in the 37T7T 0 mode, which belongs to the vector current, are contaminated by 
decays not belonging to the assigned currents (e.g. urn —> 37t), and thus need to be corrected. 
Since ~ 71 % of the 37r27r° mode consist of lotitt 0 decays, this channel is used for the unnr° 
corrections. Corrections for the other u and g modes are made with the Monte Carlo. Decay 
modes which are not reconstructed from the data have also to be included in the total vector 
and axial-vector spectra. Their distributions are taken from the Monte Carlo. The 2K7r modes 
contribute to both classes to an unknown amount. A weight of (50 ± 50) % is used for both 
currents and a correlation of —100 % between the vector and the axial-vector weights is assumed. 
The errors assigned to Monte Carlo spectra are taken to be ±100 %, in order to take a possible 
mismodelling of the Monte Carlo into account. 

The moments Ry/ A are given in table |3j. The errors on the moments are subdivided into 
statistical uncertainties due to the data statistics, the uncertainties comming from the branching 
ratio errors, and systematic uncertainties induced by the limited Monte Carlo statistics and the 
variations of the energy scale, the momentum scale, and the regularization parameter (table f|) . 
Correlations between the moments are given in table |5|. 
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moment kl 

-nkl 

ity 

total error 

TDkl 

k a 

total error 

00 

1.764 

±0.016 

1.720 

±0.017 

10 

1.264 

±0.012 

1.240 

±0.013 

11 

0.2980 

±0.0034 

0.2510 

±0.0032 

12 

0.0942 

±0.0019 

0.1090 

±0.0019 

13 

0.0403 

±0.0016 

0.0518 

±0.0013 


Table 3: The measured moments Ry/ A , for kl = 00,10,11,12,13. The errors shown represent 
statistical and systematic uncertainties. 


kl 

data stat. 

branching ratios 

MG stat. 

systematic errors 

E scale p scale regularization 


00 

- 

±0.016 

- 

- 

- 

- 


10 

±0.005 

±0.010 

±0.004 

±0.004 

±0.001 

±0.001 

V 

11 

±0.0012 

±0.0031 

±0.0008 

±0.0005 

±0.0002 

±0.0000 


12 

±0.0006 

±0.0016 

±0.0004 

±0.0007 

±0.0000 

±0.0001 


13 

±0.0008 

±0.0011 

±0.0005 

±0.0005 

±0.0001 

±0.0001 


00 

- 

±0.017 

- 

- 

- 

- 


10 

±0.004 

±0.012 

±0.002 

±0.002 

±0.002 

±0.003 

A 

11 

±0.0010 

±0.0029 

±0.0007 

±0.0003 

±0.0002 

±0.0004 


12 

±0.0008 

±0.0015 

±0.0005 

±0.0006 

±0.0001 

±0.0005 


13 

±0.0007 

±0.0008 

±0.0004 

±0.0004 

±0.0002 

±0.0003 


Table 4: Statistical and systematic uncertainties of the measured moments. The upper (lowerj 
portion of the table contains the result for the vector (axial-vector) current. 


8.2 Spectral functions 

The vector and axial-vector spectral functions are given by inverting equation (jl|): 


v/a(s ) 


2vr Im77yy A (s) 


m. 


6Sew| Vudl 



x E 


B[t —> hy/ A +-) Wy/A diVy/A 
B{t —> eu e u T ) Ay /a ds 


( 22 ) 


where the sum is performed over hadronic final states /ry/A with angular momentum J = 1. 

The spectral functions (and their correlations) are shown in figure [5| together with the flat 
naive parton model prediction u na iv e (s) = anaive(s) = 1/2 and the prediction of perturbative 
QCD (massless) for a s [m |) = 0.122 which increases the naive prediction by ~ 10 %. As a result 
of the regularized unfolding, the bin-to-bin correlations are of the order of +80% (—50%) for 
bin distances of 0.1 GeV 2 (~ 1 GeV 2 ). The correlation between vector and axial-vector spectral 
function varies from —60% to +60%. Figure || shows the difference and the sum of the two 
measured spectral functions. The function t;(s) — a(s) should vanish in the limit of perturbative, 
massless QCD. The deviation from this prediction, e.g. due to the p and ai resonances, indicates 
the large sensitivity of this distribution to non-perturbative effects. The QCD prediction for 
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-100 % 




Figure 5: The vector and axial-vector spectral functions. Shown are the sums of all contributing 
channels as data points (upper two plots). Some exclusive contributions are shown as shaded 
areas. The nawe parton model prediction is shown as dashed line, while the solid line depicts the 
perturbative, massless QCD prediction for a s (m f) = 0.122. The error bars include statistical 
and systematic uncertainties. The pion pole is subtracted from the axial-vector spectrum. The 
lower plot shows the correlations of the two spectral functions in continuous gray-levels from 
white to black which correspond to the correlations in percent from —100% to +100%. The 
contour lines are drawn in equidistant steps of 20%. 
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-7 
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12 
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0 
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13 
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8 

3 
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-26 


Table 5: Correlations between the measured moments Ry/ A in percent. The left (right) table 
gives the correlations between the moments of the vector (axial-vetor) current; the table in the 
middle shows the correlations between the moments of different currents. 


v(s) + a(s) which is ~ 10 % above the naive expectation v(s) + a(s) = 1 as in figure |] gives a 
reasonable description of the region s > lGeV 2 . The structure due to the narrow resonances 
in the region below s ~ 1 GeV 2 is however not described by perturbative QCD. 


c« 
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V 

2 

1.5 
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tB 

+ 

'x 2.5 


• OPAL 

— perturbative QCD (massless) 

— na'ive parton model 



s (GeV-) 


0 0.5 1 1.5 2 2.5 3 

s (GeV 2 ) 


Figure 6: The difference (sum) of the spectral functions v(s) — a(s) (v(s) + a(s)). The dashed 
line is the naive parton model expectation and the solid lines depict the prediction of massless, 
perturbative QCD as in figure^. For v(s) — a(s) both predictions are identically zero. 


9 Measurement of the strong coupling ct s 

Since the perturbative expansions for vector and axial-vector currents are identical while the 
non-perturbative parts have opposite sign but the same order of magnitude for both currents, 
two different fits are used for the extraction of ct s and the power corrections, respectively. The 
sum of vector and axial-vector moments is most sensitive to perturbative QCD and is used for 
the measurement of a s (fit 1) while the separate moments of both currents are used to obtain 
the power corrections (fit 2). In addition to the moments listed in table |3] it is possible to 
include the measurements of the r lifetime r T and the branching ratio = B{r —> in 
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fit 1 since each of them can be used to predict the total hadronic decay rate of the r lepton: 


R t (t t ) 

= li-1-V 

Fr F 

e ' t L e 

r„ i r„ 

(23) 

Rt(B,) 

M i M 

r e B,, r e 

(24) 


Both equations assume lepton universality so that the following equation holds: 


B 





(25) 


with if = 0.9726 

1 e 


Toll and r e = 4.0329 • 10~ 13 GeV [fTOl, [48f . The non-strange decay rate of the 
r lepton is then obtained by subtracting R T S = B s /B e = 0.150 ± 0.008 [O] from the weighted 
average R T (B fi ,r T ) of R, t (t t ) and R T {B IJ ^. In principle the electron branching ratio B e could 
also be used to determine R~ but this has a 96 % correlation with R~ from the hadronic modes 


due to the correlations of the constrained branching ratios in |TT 


Using the world average t t = (291.0 ± 1.5) fs and the fitted value B M = 0.1735 ± 0.0010 [|ll 
one gets: 

Rr(B^ T r ) - Rr,s = 3.485 ± 0.023. (26) 

From the vector and axial-vector decay rates in table one gets the following value: 


R. tN + R r \ = 3.484 ± 0.024. 


(27) 


In the hrst fit, four parameters are used to describe the five moments, leaving one degree 
of freedom for the fit: the strong coupling cx s (m^), the gluon condensate (—GG) and the 
dimension 6 and 8 contributions to the kl = 00 moments hy +A , ^v+a- The second fit requires 
six parameters to predict ten moments (four degrees of freedom): a s (m^), GG) and the 
power corrections <5y, <5y, 5^ and 5\. The power corrections from the two fits can be compared 
via the following relation: 

5v +a = \ 0° + 5?) • (28) 

Further inputs for both fits are the quark masses for the three light quarks m Uj d, s 

m u = (8.7 ± 1.5) MeV, m d = (15.4 ± 1.5) MeV, m s = (270 ± 30) MeV, (29) 

and the quark condensates ('0'0)u,d,s = — Fuds; with 

fiu = f id = (189 ± 7) MeV, /i s = (160 ± 10) MeV. (30) 

The values are taken from [Q . The error matrix of the moments is calculated from the experi¬ 
mental errors on the moments and their correlations (tables § and ||) and the theoretical error 
matrix calculated from the errors on the quark-masses and quark-condensates. The results 
from the fit to the sum of vector and axial-vector moments is given in table [|. The quoted 
errors are subdivided into a statistical error due to the data statistics, the uncertainty induced 
by the errors on the branching ratios, an experimental systematic error from the Monte Carlo 
statistics and the unfolding procedure, and a theoretical error including the uncertainties on 
quark masses, the variation of the O(cXg) coefficient, the renormalization scheme dependence, 
and the renormalization scale uncertainty. The strong coupling is most sensitive to the kl = 00 
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theory 

observable 

value 

data 

contributing errors 
B syst. 

theo. 

A 2 /d.o.f. 


a s (ml) 

0.348 

±0.002 

±0.009 

±0.002 

±0.019 


CIPT 

(fGG)/GeV 4 

-0.003 

±0.007 

±0.007 

±0.006 

±0.005 

0.16/1 


^V+A 

0.0012 

±0.0034 

±0.0033 

±0.0029 

±0.0006 



4 + a 

-0.0010 

±0.0024 

±0.0016 

±0.0015 

±0.0003 



a s (m 2 ) 

0.324 

±0.001 

±0.006 

±0.002 

±0.013 



(^GG)/GeV 4 

0.014 

±0.007 

±0.006 

±0.005 

±0.013 


FOPT 







0.17/1 


^V+A 

0.0028 

±0.0034 

±0.0034 

±0.0030 

±0.0068 



^V+A 

-0.0015 

±0.0024 

±0.0016 

±0.0014 

±0.0019 



a s (ml) 

0.306 

±0.001 

±0.005 

±0.001 

±0.011 



(^GG)/GeV 4 

-0.002 

±0.007 

±0.007 

±0.005 

±0.002 


RCPT 







0.07/1 


^V+A 

-0.0047 

±0.0036 

±0.0040 

±0.0032 

±0.0011 



^V+A 

-0.0001 

±0.0024 

±0.0017 

±0.0015 

±0.0003 



Table 6: The result for a s (rnf) and the non-perturbative parameters from the fit to the sum 
of vector and axial-vector moments. Shown are the values for the three different descriptions 
of the perturbative part of the moments (see text). The given errors correspond to the data 
statistics, the uncertainty due to the errors on the branching ratios B, a systematic error from 
the Monte Carlo statistics, the energy scale, the momentum scale, and the unfolding, and a 
total theoretical uncertainty. 


moment, and therefore the dominant contribution to the experimental uncertainty on a s comes 
from the uncertainties on the branching ratios. 

All three theories lead to similar y 2 values (see table |6j) but the spread in the fitted values 
for a s (m 2 ) exceeds the total uncertainties by a factor of two. A similar spread of the values 

im 


for aJmz) from the three models has also been observed in 


where RCPT has led to the 


lowest value and CIPT to the largest value in agreement with our results (table [|) . 

The differences in the statistical and systematic errors on ct s are induced by the scaling of the 
relative error with ct s and thus are compatible for the three fits. The theoretical uncertainties 
should also obey this scaling behavior: here the fits for FOPT and CIPT only include the 
uncertainty on the unknown Jt 4 coefficient and hence cannot be compared to the RCPT result. 
Furthermore the uncertainty due to the variation of the renormalization scheme vanishes for 
RCPT. The impact on a s from the various theoretical error sources is listed in table jTj. The 
given errors correspond to the spread of the fitted values of a s in fit 1 due to the unknown 
O (ccf) dependence K 4 = 25 ± 50, the choice of renormalization scale 0.4 < p 2 /m 2 < 2.0, 
the variation of the renormalization scheme parameterized with the third coefficient of the 
/3-function 0.0 < l3f’ s //3^ is < 2.0, and the evolution of a s (m 2 ) to the Z°-mass scale. 

Although the total theoretical uncertainties on a s are compatible for all three theories 
there is a major difference between FOPT and the two other models: the FOPT fit leads to 


6 R T 'S is subtracted from R T , since the induced dependency on the mass of the strange quark would lead to 
a larger uncertainty in the fits if R T would be used instead. 
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error source 

CIPT 

Aa s (m 2 ) 

FOPT 

RCPT 

CIPT 

Aa s (m|) 

FOPT 

RCPT 

-25 < K a < 75 

±0.012 

±0.006 

- 

±0.0013 

±0.0007 

- 

0.4 < p 2 fm 2 < 2.0 

±0.006 

±0.009 

±0.011 

±0.0005 

±0.0009 

±0.0015 

0.0 < /3^ s //3p < 2.0 

±0.015 

±0.009 

±0.000 

±0.0009 

±0.0005 

±0.0005 

evolution 

- 

- 

- 

±0.0003 

±0.0003 

±0.0003 


Table 7: The theoretical uncertainties on the strong coupling constant. The errors correspond 
to the full spread of the fitted a s values in fit 1 due to the variation of the parameters listed in 
the first column. 


a significant larger dependency of the non-perturbative parameters (-GG) and (fy+A on the 
theoretical uncertainties than CIPT and RCPT. The dominant effect comes from the variation 
of the renormalization scale p 2 . The statistical and systematic uncertainties on the power 
corrections are very similar for all three theories, agreeing with expectation. Figure [7] shows 
a comparison of i? r v(s 0 ) + R T! a(sq) as predicted from the three theories using the fit results 
at So = m 2 with the data. The Contour Improved prediction is consistent with the data from 
the T-mass scale down to so ~ 1 GeV 2 while FOPT and RCPT tend to predict too large values 
below s 0 ~ 2 GeV 2 . 


9.1 Evolution of cx s from m T to mz 


The value of the strong coupling at the mass scale of the r lepton can be evolved up to the 
mass scale of the Z°. This is done by solving the four-loop /3-function given by equation (7|) 
numerically in small steps from rn 2 to m| applying a three-loop matching condition 


m 


flavor thresholds for m c (m c ) = (1.30 ± 0.06) GeV and m h (m h ) = (4.13 ± 0.06) GeV fl50 
evolution procedure induces an additional error of ±0.0003 [ 
mass. Using the CIPT result for a s {m 2 ) and mz = 91.187GeV the following value is obtained: 


at the 

. The 

on the strong coupling at the Z° 


a s (m|) 

= 0.1219 ± 0.0010e X p ± 0.0017 theo ± 0.0003 evol . 

(31) 

The FOPT fit gives 



a s (m 2 z ) 

= 0.1191 ± 0.0008 exp ± 0.0013 theo ± 0.0003 evol . 

(32) 

Finally RCPT gives: 



a s (m 2 z ) 

= 0.1169 ± 0.0007 exp ± 0.0015 theo ± 0.0003 evol . 

(33) 


The different contributions to the theoretical uncertainties are listed in table [?]. The results are 
in good agreement with the value obtained from fits to combined electroweak measurements at 


LEP and SLD 15TI 


a 8 (m|) = 0.120 ±0.003. 


(34) 
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Figure 7: The non-strange hadronic decay rate of the r lepton R Tt v(s 0 ) + R t ,a(so ) versus the 
upper integration limit s 0 . The points in the upper plot denote OPAL data; the dashed, dashed- 
dotted and dotted curves represent the theoretical predictions of the three theories with the results 
from the fit to the moments at s 0 = mf used as input. The lower plot shows the three theories 
normalized to the data. The three sets of dashed, dashed-dotted and dotted curves indicate 
central values and total experimental errors for each of the three theories. The errors on the 
data are shown as solid curves. 
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contributing errors 



theory 

observable 

value 

data 

B 

syst. 

theo. 

X 2 / d.o.f. 


ot s (ml) 

0.347 

±0.001 

±0.012 

±0.002 

±0.019 



{fGG)/GeV 4 

0.001 

±0.003 

±0.006 

±0.003 

±0.004 



4 

0.0256 

±0.0017 

±0.0024 

±0.0017 

±0.0006 


CIPT 

4 

-0.0080 

±0.0010 

±0.0007 

±0.0005 

±0.0002 

0.63/4 


% 

-0.0197 

±0.0016 

±0.0022 

±0.0019 

±0.0010 



si 

0.0041 

±0.0012 

±0.0013 

±0.0008 

±0.0002 



(ml) 

0.323 

±0.001 

±0.008 

±0.002 

±0.014 



(fGG)/GeV 4 

0.017 

±0.003 

±0.004 

±0.003 

±0.010 



5 6 V 

0.0271 

±0.0017 

±0.0025 

±0.0018 

±0.0056 


FOPT 

S 8 y 

-0.0085 

±0.0010 

±0.0007 

±0.0005 

±0.0012 

0.62/4 


si 

-0.0183 

±0.0016 

±0.0023 

±0.0019 

±0.0052 



si 

0.0036 

±0.0011 

±0.0012 

±0.0008 

±0.0011 



a s (m 2 r ) 

0.305 

±0.001 

±0.007 

±0.001 

±0.011 



{fGG)/GeV 4 

0.002 

±0.003 

±0.005 

±0.003 

±0.001 



si 

0.0202 

±0.0018 

±0.0033 

±0.0018 

±0.0009 


RCPT 

si 

-0.0075 

±0.0010 

±0.0008 

±0.0005 

±0.0002 

0.61/4 


Si 

-0.0252 

±0.0017 

±0.0032 

±0.0020 

±0.0006 



si 

0.0047 

±0.0012 

±0.0013 

±0.0008 

±0.0001 



Table 8: The fit result for a s and the power corrections from the combined fit to vector and 
axial-vector moments. The given errors correspond to the data statistics, the uncertainty due 
to the errors on the branching ratios B, a systematic error from the Monte Carlo statistics, the 
energy scale, the momentum scale, and the unfolding, and a total theoretical uncertainty. 


10 Measurement of dimension 6 and 8 operators 

The results from fit 2 where the separate moments of the vector current and axial-vector current 
are used are given in table || In contrast to a s where the error is dominated by the theoretical 
uncertainties, the power corrections are almost independent of the theoretical uncertainties for 
CIPT and RCPT. As mentioned in section [|, this is not the case for the FOPT fit which leads 
to theoretical errors of the order of (or even larger than) the experimental errors. Due to the 
correlated unfolding of vector and axial-vector spectra a strong positive correlation between 
the power corrections of the vector and axial vector current of the same dimension is observed. 
The power corrections of different dimension but for the same current are anti-correlated. All 
correlations of the fit parameters for CIPT are summarized in table |3. The fitted values of 
the strong coupling constant in both fits are in excellent agreement for all three models. The 
experimental error on a s from this fit is larger than in fit 1 as the additional information from 
the r lifetime and the branching ratio B{r —> /iz/ M z/ T ) is omitted. Using equation (|28|), the 
separate and total power corrections are also in good agreement for all three models. As in fit 
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Table 9: Correlations between the QCD parameters from the fit to the moments of the vector 
and axial-vector current in percent. The given numbers are taken from the CIPT fit result. 


1 all three theories give similar y 2 values in the fit to the exclusive moments. The theoretical 
uncertainties behave similarly in fit 1 and fit 2. The sum of all power corrections <5 n on-pert,v/A 
and <5non-pert,v+A to R r ,v / a and R r y ± R T) \ including the dimension 2 quark-mass correction 
and the dimension 4 correction obtained from the fitted gluon condensate are: 


0.0172 ±0.0026 

CIPT 


<W P ert,v = 0.0187 ± 0.0054 

FOPT 

(35) 

0.0124 ±0.0033 

RCPT, 


-0.0219 ±0.0026 

CIPT 


5non-pert,A = —0.0204 ± 0.0050 

FOPT 

(36) 

-0.0266 ±0.0032 

RCPT, 


-0.0024 ±0.0025 

CIPT 


<Wpert,V+A = -0.0009 ±0.0051 

FOPT 

(37) 

-0.0071 ±0.0031 

RCPT, 



where the errors include experimental and theoretical uncertainties. Thus all three theories 
lead to non-perturbative corrections to R r y (R t ,a) of the order 1.6% (—2.3%), while a large 
cancellation of both contributions leads to a total non-perturbative correction to R r y ± R T ,A 
which is compatible with zero and therefore allows a precise measurement of the strong coupling 
constant in fit 1. The numbers in table || can be compared to the estimates given in 0: 


(—GG)/GeV 4 = 

7T 



0.02 ± 0 . 01 , 

0.024 ±0.013, 
-0.038 ±0.020, 
- 0 . 0001 . 


(38) 


Only the power corrections of dimension 8 seem to be underestimated, while the other estimates 
are in good agreement with the measured values. Figure |5] shows the two power corrections of 
dimension 6 (CIPT) together with the theoretical prediction given in |4[]. 
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Figure 8: The power corrections of dimension 6 to R t ,v/a■ Shown are the one and two stan¬ 
dard deviation contours of the correlated result for the vector and axial-vector current (CIPT) 
including experimental and theoretical uncertainties. The solid line indicates the theoretical 
prediction given in ||]. 


11 Test of the ‘running’ of ct s 

The fit to the sum of vector and axial-vector moments (fit 1) can be extended to lower values 
of So, thus giving a correlated measurement of the strong coupling at different scales. Four 
equidistant values for So between 1.3 GeV 2 and m 2 are used. In addition to the five moments at 
so = m 2 the integrated differential decay rate i?°° v (so) + i?°° A (so) for each additional so value 
is included in the fit (see figure |?i) . 

For the extraction of the ‘running’ of a s the number of fit parameters is increased to include 
the strong coupling ct s (so) for each so value below m 2 . The result can be examined with the 
four-loop /3-function. This is shown in figure [|. where the /3-function has been refitted for all 
three sets of a s values. The values at so = 1.3 GeV 2 were not included in the fit. A comparison 
of these values with the predicted ‘running’ shows good agreement in case of CIPT, while a 
weaker ‘running’ as predicted by the /3-function is preferred by the FOPT and the RCPT values. 


Figures [7j and can be regarded as tests of the validity of the OPE for so values below 
m 2 . It has been questioned if the definition of R T (so) is still valid in this region ||, since the 
endpoint s = s 0 is no longer suppressed by the (1 — s/m 2 ) 2 term in front of the spectral function 
(see equation ([!])). By defining the hadronic decay rate for a hypothetical r' with a mass of 
m T / = yGo and inserting m T > for m T in equation (jT) one gets ||: 


so 

Rr',v/ a(sq) = 127rS'Ew|Kid| 2 y 


ds 

s o 


1 - -)' 
So/ 


1 + 2 - 


s 0 


Im-ff+l 


s) + ImlJyj A (s) 


(39) 


obeying the same quadratic suppression of the endpoint on the real s-axis as R Tt y/\(ml). 
Figure [B] shows the sum R r iy(so) + R T > t A(so) versus the upper integration limit Sq. The error 








S 0 (GeV 2 ) 


Figure 9: The ‘running ’ of the strong coupling. The three sets of a s values are shown as 
data points. The error bars include statistical and systematic uncertainties. The dashed curves 
represent the predictions of the 4-loop /3-function obtained from fits to the three sets of a s values 
not including the values ct s (l .3 GeV 2 ). The solid lines depict the errors from the fits. 

band for CIPT in the lower plot shows that the uncertainties increase below so — 1.5 GeV 2 
compared to the error in the lower plot of figure |7|. While the error on R t (sq) is dominated 
by the uncertainty of the perturbative expansion, the error on R t i(sq) originates mainly from 
its dependency on the non-perturbative parts. In contrast to R T where these power corrections 
stay constant for all s 0 (see equation m they increase with powers of 1 /s 0 as s 0 decreases 
in the case of R r >. As the errors are large for small values of so little can be said about this 
region. 


12 QCD sum rules 


Weighted integrals over the difference of the two measured spectral functions shown in figure joj 
can be compared to the chiral predictions of QCD sum rules: 


h(s 0 ) = ^2 f ds ( u ( s ) - a ( s )) = 

o 

1 S ° 

h(s 0 ) = ^Jdss (v(s) - a(s)) = 0, 
o 

1 7 ds 

J — (v(s) - a(s)) = ft ^ - F A , 


hysf) — 


(40) 

(41) 


c2 (rl) 


47T 2 


29 


3 


(42) 








Figure 10: The non-strange hadronic decay rate of a hypothetical r' lepton with rnf, = so versus 
the upper integration limit s 0 . The points in the upper plot denote OPAL data; the dashed, 
dashed-dotted and dotted curves represent the theoretical predictions of the three theories with 
the results from the fit to the moments at s 0 = mf used as input. The lower plot shows the 
three theories normalized to the data. The three dashed curves indicate central values and 
total experimental errors for CIPT. The dashed-dotted and dotted curves show central values 
for FOPT and RCPT. The errors for FOPT and RCPT are similar to the CIPT errors and 
omitted from the plot. The errors on the data are shown as solid lines. 
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S 0 (GeV-) 


Figure 11: QCD sum rules. Equations are shown in the plots (a)-(d). Shown are 

the integrals versus the upper integration limit as shaded bands. The chiral prediction is given 
by the lines (±1 a when two lines are present). 


SO 


h(s 0 ) = ^ jdss ln^ O(s) - a(s)) = 
o 


(«± “ m l°) • 


(43) 


Here the right hand side of each equation is understood to be the chiral prediction in the limit 
oo. Equation 


So —* oo. equation ( ]4U| ) is known as the first Weinberg sum rule [52], assuming that the only 
scalar contribution is given by the pion pole which is related to the pion decay constant f n = 
(92.4±0.26) MeV [|llj]. The second Weinberg sum rule [|52| is given in equation (f4T|) . The Das- 
Mathur-Okubo (DMO) sum rule [ fa3| is given by equation ( |42| ) . Its asymptotic prediction is a 
function of the pion decay constant f n , the mean square of the pion charge radius (r^) = (0.439± 
0.008) fm 2 |54j and the axial-vector form factor of the pion F A = 0.0058 ± 0.0008 [|lTI[], and 
equation ( |45| ) gives the electromagnetic mass difference of pions [p3| . Note that equation ( |45| ) 
does not depend on the cut-off value A by virtue of the second Weinberg sum rule. 


7 Our definitions of F A and differ by a factor of 1/2 from those given in 
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The saturation of these four sum rules is tested taking into account the full correlations 
between the measured spectral functions. The plots of figure pd] show the measured values 
of the integrals / 1-/4 as error bands including all experimental uncertainties versus the upper 
integration limit. The asymptotic predictions are drawn as thin lines denoting the present ±1 a 
ranges. 

All four sum rules appear to be saturated at the r-mass scale within their errors. However, 
due to the small phase space near the r-rnass which appears in the denominator of the spec¬ 
tral functions these errors are very large except for the DMO sum rule where the 1/s factor 
suppresses the high energy tail. The value of equation ([12]) at s 0 = rn 2 is: 

h(m 2 T ) = (26.3 ± 1.8) • 10" 3 , (44) 

where the error covers all experimental uncertainties. 


12.1 Pion polarizability 

Assuming that the DMO sum rule shown in figure 0(c) is already saturated at the r-mass 
scale, its value can be used to predict the electric polarizability of the charged pion as proposed 

in Q: 

_ _a_ ({rp h\ , . 

E m,* ( 3 f}) ' ( ) 

Using the result from the previous section for the DMO sum rule (equation (f|4|)) one gets: 

a E = (2.71 ± 0.88) • 10 " 4 fm 3 , (46) 

which is in good agreement with the value a E = (2.64 ± 0.36) • 10 _ 4 fm 3 , derived in |56 . 


13 Summary 

Measurements of the spectral functions of the vector current and the axial-vector current and 
their applications in QCD have been presented. Within the framework of the Operator Prod¬ 
uct Expansion, a simultaneous determination of the strong coupling constant a s and non- 
perturbative correction terms has been performed. The sum of R r y and R T) a was found to 
involve a large cancellation of the non-perturbative terms and thus has been used together with 
the r lifetime and the branching ratio B(r —> /i/z /t /z T ) to give a precise measurement of the 
strong coupling constant. CIPT has led to the value 

a s (m 2 r ) = 0.348 ± 0.009 exp ± 0.019 theo , 


at the r-mass scale and 


a s (m 2 z ) = 0.1219 ± 0.0010 exp ± 0.0017 theo 

at the Z°-mass scale, where the first error stems from the experimental uncertainties and the 
second error originates from the theoretical uncertainties. The values obtained for cc s (m|) using 
FOPT or RCPT are 2.3% and 4.1% smaller, respectively. 

The total amount of non-perturbative corrections to R, r y (7? t ,a) was found to be (1.6±0.4) % 
((—2.3 ± 0.4)%), while the correction on the sum of R T y and R T) a due to non-perturbative 
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QCD is found to be only (—0.3±0.4) %. Here the errors include all experimental and theoretical 
uncertainties. 

Assuming the validity of the Operator Product Expansion for energy scales below the r 
mass a test of the ‘running’ of the strong coupling between so — 1.3 GeV 2 and so = m 2 has 
been performed. A good agreement between the predicted ‘running’ from the 4-loop /1-function 
and the fitted a s values has been observed for CIPT. 

The saturation of QCD sum rules at the r-rnass scale has been tested, yielding a measure¬ 
ment of the pion polarizability of «e = (2.71 ± 0.88) • 10~ 4 fm 3 as determined from the DMO 
sum rule. 
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